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S-wave charmed mesons in lattice NRQCD 
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Heavy-light mesons can be studied using the 1/M expansion of NRQCD, provided the heavy quark mass is 
sufficiently large. Calculations of the S-wave charmed meson masses from a classically and tadpole-improved action 
are presented. A comparison of 0(1/M), 0(1/M 2 ) and 0(1/M 3 ) results allows convergence of the expansion to 
be discussed. It is shown that the form of discretized heavy quark propagation must be chosen carefully. 
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1. INTRODUCTION 

Nonperturbative strong dynamics is typified by 
the mass scale Aqcd- Interactions involving a 
very heavy quark of mass M can be studied sys- 
tematically by expanding in Aqcd/M. Upon 
truncation of the expansion at some order, the re- 
sulting effective theory is not renormalizable and 
requires a momentum cutoff. 

If the required regularization is performed via a 
space-time lattice, then the cutoff is proportional 
to the inverse lattice spacing 1/a. A useful effec- 
tive theory must satisfy 

Aqcd <C 1/a < M, (1) 

so the cutoff is large enough to include the bulk of 
the nonperturbative dynamics in the low-energy 
effective theory, but small enough that the trun- 
cation of the expansion remains sensible. 

Is the charm quark heavy enough for a use- 
ful lattice effective theory? The present work 
addresses this question through a study of 
the masses of S-waves charmed mesons using 
quenched lattice NRQCD. [0 Calculations are 
performed at two lattice spacings near 0.22 fm 
and 0.26 fm, and in each case results are given 
separately at 0(1/M), 0{1/M 2 ) and 0(1/M 3 ) 
in the effective theory. This work is an extension 
of results that have been reported previously 
and further details can be found in that paper. 
Other authors have considered NRQCD up to 
0(1/M 2 )§. 



2. ACTION 

The lattice action has three terms: gauge ac- 
tion, light quark action and heavy quark ac- 
tion. The entire action is classically and tadpole- 
improved with the tadpole factor defined by 
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The gauge action includes a sum over 1x2 rect- 
angular plaquettes as well as 1 x 1 elementary pla- 
quettes. For light fermions, the Sheikholeslami- 
Wohlert action^ is used with the clover coef- 
ficient set to its tadpole-improved value. The 
heavy quark action is NRQCD. 

A discretization of the NRQCD action leads to 
the following Green's function propagation Q: 
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2n J Uq V 2n , 
x(l-a5H)G T , t>0, (4) 

where "n" should be chosen to stabilize the nu- 
merics, and the Hamiltonian is 



H 
SH 

Ho 
SHU 



H a + 5H, 
SHU + 6H& 



(5) 

0(1/Af 4 ), (6) 



2M 
c 4 



L7 4 2M 



B 



C5- 



t 2 A' 4 ' 
24M : 



(7) 

(8) 



2 



SH^ 2 



c 2 w , : - - z , a(A( 2 )) 2 

-f^-( AXE - EXA )' ( 9 ) 



-Ci- 



(A 



(2)n,2 



C7 g 

8Af 3 U$ Ml 3 



C/ 8 8Af 3 



er • (E x E + B x B) 



™ <* (E» + B')-c 1I ° 2 ' a(2) > 3 



C/ 8 8Af 3 



192n 2 M 3 

+purely quantum effects, (10) 
F 4l , (11) 



1 



(12) 



A tilde denotes removal of the leading discretiza- 
tion errors. Classically, the coefficients c$ are all 
unity, and their nonclassical corrections will not 
be discussed in this work. 

It should be noted that the separation of Hq 
and SH in Eq. (^) is not unique. The heavy quark 
propagation of Eq. (^) uses a simple linear ap- 
proximation to the true exponential dependence 
on SH , while using a better-than- linear approxi- 
mation for Hq. The present work will report on 
a generalization of this choice. 

3. RESULTS 

All data presented here correspond to a 
charmed meson with a light quark mass that is 
roughly twice the strange quark mass. More 
extensive results, including bottom mesons, can 
be found in Ref. ||. The data sample in- 
cludes 400(300) gauge field configurations at /3 = 
6.8(7.0) corresponding to a « 0.26fm(0.22fm). 
All plots include bootstrap errors from 1000 en- 
sembles. 

Fig. [l] shows the simulation energy (as read 
from the plateau of an effective mass plot) of 
the 1 S'o meson. Notice that at both lattice spac- 
ings the 0(1/M 3 ) contribution is twice as large 
as the 0(1/M 2 ) contribution, and that this large 
effect is dominated by the term containing cio 
in the Hamiltonian, Eq. (|To|). This term is 
unique because it is the only term up to 0(1/M 3 ) 
which contains a nonzero vacuum expectation 



value. The vacuum value can be calculated from 
our gauge field configurations, and as shown in 
Fig. |l| the simulation energy displays a very pleas- 
ing 1/M expansion after removal of the vacuum 
value. 

The mass difference between 3 Si and 1 So 
mesons is shown in Fig. |[ When the vacuum ex- 
pectation value is not removed from the Hamilto- 
nian, the 0(1/M 3 ) contribution is twice as large 
as the 0(1/M 2 ) contribution in magnitude, due 
to large effects from the terms containing C7 and 
cio . Fig. H indicates that the cio effect is entirely 
due to the vacuum expectation value. 

This very substantial dependence of the spin 
splitting on the vacuum value should be disturb- 
ing, since the c\q term is spin- independent. Ap- 
parently the vacuum value is so large that it 
destabilizes heavy quark propagation and thereby 
introduces a spurious effect unless the vacuum 
value is removed from the Hamiltonian. To sup- 
port this claim, we have redone the calculation af- 
ter subtracting the cio term from SH and adding 
it to Hq in Eq. (Q). When this is done, the trian- 
gles of Fig. Hare reproduced regardless of whether 
the vacuum value is subtracted from the Hamil- 
tonian. This is the physically-expected result. 

Still, Fig. H contains a large 0(1/M 3 ) con- 
tribution dominated by the term containing c-j. 
Although this term does not contain a nonzero 
vacuum expectation value, one wonders if the 
heavy quark propagation might be unstable for 
this term as well, unless a better-than-linear ap- 
proximation is used for the C7 term in Eq. (|j). 

Fig. HI shows the effect of subtracting all of SH 
from its present location in Eq. (^) and putting 
the full Hamiltonian in place of Hq. That is, 
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and we choose Go = fe.o- To maintain classical 
improvement, the following (3(1/M 3 ) term must 
be added to the Hamiltonian: 

SH ncw = -± { (h + 5HW) , SH^ } . (14) 

Simulations with n = 5 and n = 7 are indistin- 
guishable. At both lattice spacings, the large C7 
effect is found to be robust, and the contribution 
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from eg tends to increase. No discussion of the 
spin-independent terms containing c\ and en is 
presented here. 

4. DISCUSSION 

Instabilities can arise in the NRQCD expan- 
sion for charmed mesons due to the presence of a 
large vacuum expectation value. A better-than- 
lincar approximation to heavy quark propagation 
is valuable for ensuring stability. 

Substantial effects on spin splitting were found 
from 0(1/M 3 ) terms in the action. However fur- 
ther study, for example, of alternative definitions 
for the tadpole factor or of perturbative improve- 
ment for the NRQCD coefficients is needed be- 
fore one can reach a definitive conclusion about 
the convergence of the NRQCD expansion for 
charmed mesons. 
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Figure 1. The simulation energy of a ground 
state charmed meson at rest, up to 0(1/M k ). 
Solid(open) symbols denote data at (3 = 6.8(7.0). 
Triangles are produced by subtracting the vac- 
uum expectation value from the cio term. To the 
right of the vertical line, the effect of adding each 
0(1/M 3 ) term to the 0(1 /M 2 ) Hamiltonian is 
shown individually. 
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Figure 2. The spin splitting of S-wave charmed 
mesons, with notation as in Fig. Ill 
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Figure 3. The data (circles) of Fig. |, which use 
Eq. (^) , are compared to data (squares) obtained 
from Eq. (|l3[). The vacuum expectation value is 
not removed. 



